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6, Theorie van zonnevlammen 


M. Kuperus 


We zullen in dit college uitsluitend die theorieën bespreken, die er 
van uitgaan, dat de energie die bij een vlam vrijkomt energie is die in 
magnetische velden is opgeslagen. Het probleem van de zonnevlam valt dan 


uiteen in vier deelproblemen, nl,: 


1. de opslag van energie in magnetische velden; 


[na 
e 


de energievrijmaking door een magnetohydrodynamische instabiliteit; 


de reconnectie van magnetische velden en de vorming van een veldcon- 


fu» 


figuratie met een lagere energie; 


de versnelling van deeltjes die tijdens de processen in 2) en 3) op- 


je 


treedt. 


Een potentiaalveld B e= Vb is de toestand met minimale energie als 
de potentiaalverdeling op de rand van een systeem is gegeven. Het is dui- 


delijk dat in een potentiaalveld geen stromen lopen omdat 
i= gVXB e= VXV e= 0, (1) 


Ieder niet-potentiaal veld bevat meer energie, Deze energie of vrije 
energie zit precies in de stromen j . In vacuum kunnen geen stromen lo- 
pen en daarom is ieder vacuum veld een potentiaal veld. 


De magnetische energie in een volume V is 
B? 
Ig > |E dv, (2) 


De enige manier om W, te vergroten is B te vergroten. Dit kan op twee 


manieren gebeuren nl. 
&. door de veldsterkte op de rand te vergroten 
(eigenlijk de flux vergroten); 


b. door zonder fluxvergroting de polariteiten anders te verdelen. 


In het eerste geval zal het duidelijk zijn dat weer een potentiaalveld 
kan worden verkregen echter met een grotere energie=inhoud. 
Dit nieuwe veld bezit evenwel geen vrije energie. De energietoevoer zit 


in het mechanisme dat B op de rand heeft vergroot. 


In een plasma gaat vergroting van de veldsterkte (of flux) altijd 
gepaard wet compressie als het plasma zeer goed geleidend is en alleen 
bewegingen die met compressie gepaard gaan kunnen zo tot een potentiaal- 
veld met hogere energie aanleiding geven. 

In het tweede geval heeft men alleen met redistributie van de pola=- 
riteiten te maken of strikter met beweging van de voetpunten van de 
krachtlijnen. We moeten goed bedenken dat in een ideaal plasma (cq = °9) 
de krachtlijnen niet slechts een mathematisch concept zijn maar een 
materiële betekenis hebben waarmee zeer goed kan worden gewerkt. In het 
geval b) is het niet zonder meer duidelijk of de energie toeneemt of 
niet, De uiteindelijke potentiële veldconfiguratie zal evenveel energie 
bevatten als de aanvankelijke, Wel zal tussen de begintoestand en de 
eindtoestand die beide een potentiaalveld zijn het veld zich in een hoge- 


re energietoestand kunnen bevinden: 





figuur | 


In bovenstaande figuur is schetsmatig aangegeven hoe de magnetische ener= 
gie verandert als men de flux verhoogt (a) en als men de flux constant 
houdt maar de voetpunten beweegt (b). De getrokken lijn geeft de energie 
in een niet-potentiaal veld en de onderbroken lijn die in een potentiaal 
veld. Er is aangenomen dat og ce zodat tenslotte de energie in het veld 
die van een potentiaal veld wordt. De tijd to is de tijd gedurende welke 
het mechanisme van energietoevoer aan het veld werkzaam is: de opslagtijd, 
terwijl de tijd tj de karakteristieke tijd is voor het afvloeien van de 


magnetische energie. Om energie in een vled te kunnen opslaan moet gelden 
tr > to. (3) 


Om echter in een korte tijd de vrije energie uit het veld te krijgen (dus 
At << to) moet zich een bijzondere situatie voordoen die we een instabi- 


liteit noemen. We gaan daar in 5 2 op in. 


di 


De toename van We wordt beschreven door (2) naar t te differentie- 
ren en gebruik te maken van de MHD-vergelijking 
B 


zE = VXxyxB+ (40) vB. (4) 


Vergelijking (2) wordt dan 


hi 


rade =| { B, (Vv x vx B) + (áro) } B. v? B} dv . (5) 


Nu geldt met behulp van de vectorrelaties: 


| crop. ven ov = - grot fe.vxp es = 


Sl- 


= — (4mo) * | [ 7. (jxB) +(axB).jl dv = 
n -[Ew, (6) 


omdat de eerste term tussen vierkante haakjes verdwijnt na integratie 


over V hetgeen men inziet door de stelling van Green toe te passen: 


[v.aar= fa.as =0 als A=0 op S. (7) 
v IS) 
De tweede term in (5) stelt de energie voor die in warmte is omgezet 
t.g.v. het eindige geleidingsvermogen. 
Op dezelfde manier kan men de eerste term in (5) m.b.v. de stelling van 
Green omzetten in 


Lü 


Ef v.tenard erm. vena « faxp.va. 
Vv 


v (8) 
Vergelijking (5) wordt nu 

og : de 
era fax» ear 5ov, (9) 

v 

waarin de eerste term de arbeid is die moet worden verricht om met een 
snelheid v tegen de Lorentzkracht j Xx B te werken, De eerste term 
stelt de energietoevoer aan het magnetische veld voor en de tweede term 
de energieafvoer. Let wel op dat (9) is afgeleid onder de veronderstel- 
ling dat j = 0 buiten een vast oppervlak. Als de stromen zich tot bui- 
ten dit oppervlak uitbreiden komt er nog een derde term bij die het ener- 
gieverlies door golven naar buiten voorstelt. Wij zullen dit effect voor- 
lopig verwaarlozen en nu enige gevallen bespreken waarbij energie-opslag 


in coronale magneetvelden mogelijk is. 


Uit vgl. (9) volgt direct dat energie=opslag zeer goed mogelijk 
wordt als het plasma goed geleidend is (zoals in de corona) en het ge- 
bied waarbinnen de stromen vloeien beperkt blijft. 

Het eerste geval is dat van het uitrekken van een bipolair magne=- 


tisch veld door radiale stroming zoals hijv. de zonnewind. 





figuur 2 


Bij het uitrekken ontstaat een stroomlaag die ken overgaan in een neutra- 
le laag. Hwr openen van de veldlijnen wordt ten naaste bij bepaald door 
de hoogte waarop de snelheid gelijk wordt aan de Alfvêén snelheid 
en A J hé 
(470) 


Deze veldconfiguratie vertoont veel overeenkomst met de magnetosfeer. 


Het tweede geval is het opwikkelen van een magnetische veldbuis door 


fotosferische bewegingen. 





figuur 3 


Hier wordt een stroom langs de veldlijnen opgewekt zodat behalve een 
axiale component B, een azimuthale component B, aanwezig is. De 
extra energie zit in het azimuthale veld. 

Het derde geval is dat van een roterende magnetische vlek in een 


omgeving van tegengestelde polariteit. 





figuur à 


Ook nu wordt het veld opgewikkeld en zal het zich ontwikkelen via een 
‚aantal evenwichtstoestanden zoals krachtvrije velden, waarop we nog terug 
komen, 

Al deze configuraties hebben model gestaan in zonnevlamtheorieën; 
het eerste model (A) in de theorie van Sturrock, het tweede model (B) in 
de theorie van Gold en Hoyle en in die van Alfvén en Carlquist en het 
derde model (C) in de theorie van Barnes en Sturrock. We zullen later nog 
laten zien dat in alle modellen tenslotte een sterke stroomlaag optreedt 
waar de dissipatie van de in het veld opgeslagen energie plaatsvindt, 

Naast de bipolaire magnetische veldconfiguraties zijn er ook theo- 
rieën die gebruik maken van quadrupool velden. Een voorbeeld hiervan is 
de theorie van Sweet, waar twee bipolaire magnetische gebieden naar el- 
kaar worden toebewogen, Wanneer die beweging snel plaatsvindt ontstaat 
een magnetische neutrale laag waarin de energie, nodig om het veld te 


verplaatsen, wordt opgeslagen. 


\ 





figuur 5 


Het verschil met de vorige theorieën is dat nu de voetpunten in bipolaire 
paren t.o.v. elkaar worden bewogen. Van energie-opslag is pas sprake als 
de beweging van het ene bipolaire gebied wordt bemoeilijkt door aanwezig- 


heid van het andere. 


6 


De theorie van Syrovatskii over dynamische dissipatie is een variant 
van de theorie van Sweet, waarbij in plaats van twee bipolaire gebieden 


gebruik wordt gemaakt van twee evenwijdige stromen (in punten A). 





figuur 6 


In punt B is een X-type neutraal punt. De veldconfiguratie in boven- 
staande figuur is het potentiaal veld van twee evenwijdige stroomdraden 
met gelijke stroomsterkte, 

Let er op dat de veldrichtingen in het neutrale punt loodrecht op elkaor 
staan. Een loodrechte stand van de veldasymptoten in B is kenmerkend 
voor een minimale energie. 

Veronderstel dat nu de stromen in A naar elkaar toe bewegen, Dat 
is niet zo moeilijk omdat gelijkgerichte stromen elkaar aantrekken, Dan 
zal bij langzame toenadering dezelfde veldconfiguratie blijven bestaan 
als in figuur 6. Bewegen we de stromen sneller dan de diffusiesnelheid 
naar elkaar toe, dan kan reconnectie niet snel genoeg plaatsvinden en 


hopen plasma en magnetisch veld zich op ín een neutrale laag. 
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figuur 7 


De neutrale laag BjB2 bevat in feite een stroomsysteem dat de toenade- 


ring van de stromen in A tegenwerkt. De invloed van het plasma is hier 


een resistieve. 

Het zal een ieder duidelijk zijn dat in stroomlagen veel energie kan 
worden opgeslagen. Toch zijn er meer bescheiden configuraties waarin ook 
energie kan worden opgeslagen, maar in minimale hoeveelheden: dit zijn de 
2,g. krachtvrije velden waarvoor geldt dat overal de Lorentzkracht 
ixB= 0. p 

In feite kan men alleen van een krachtvrij veld spreken als Pe << Br 
dus als B << 1 omdat anders aan de magnetohydrostatische vergelijking 


moet worden voldaan. 


F3 
Vp + gr (VXB) XB = 0. (10) 


Velden die aan bovenstaande vergelijking voldoen heten ook wel Lundquist 
velden. 
Als B << Ì wordt (10): 


(Vv xB)xB = 0, (11) 
waaruit volgt 

VxB = B, (2) 
waarbij ,4(r) een constante is langs de veldlijnen, maar van veldlijn tot 
veldlijn kan veranderen. a geeft de grootte aan van de stroomdichtheid 
j langs het magnetische veld. 


Splitst men het veld in een axiaal veld B, en een azimuthaal veld 


dan geldt voor B 


B» Ô 
27 
| B, rdò = If j ds . (13) 
o mr? 

Dus B, = 5 (14) 


ro° j 
8 





ne (15) 
r 


voor r > ro als j = 0 voor r > ro. 
Het azimuthale veld B, valt dus af met z buiten een stroomfilament en 
het is daarom niet mogelijk een krachtvrij filament of een krachtvrije 
fluxbuis als zelfstandige eenheid met radius ro te beschouwen en te ver- 
onderstellen dat daarbuiten geen magnetisch veld is. 

Om over een fluxbuis te kunnen spreken heeft men Öf een stroomlaag 
rond de buis nodig öf een voldoende hoge gasdruk buiten de buis, zodat 
vgl. (10) kan gelden, 


hifi Teho alen in nil nd ohne 


De ‘vrije! energie zit in het veld B, \ 


We |t ev. (16) 


Een fundamentele vraag is onder welke omstandigheden een niet potentiaal 
veld stabiel is. 

Het ligt voor de hand deze vraag eerst te stellen voor de kracht- 
vrije velden. 

Als krachtvrije velden stabiele velden zijn, is het mogelijk doer 
een serie van evenwichtstoestanden in een hogere energetische toestand 
te komen zonder dat door een instabiliteit de extra energie het systeem 
weer verlaat. Dit blijkt geen eenvoudig probieem te zijn en slechts in 
bepaalde gevallen kan hierop een antwoord worden gegeven. Er zijn enkele 
manieren om dit probleem aan te pakken. Het komt er in het algemeen op 
peer dat men een energie-integraal W opstelt en een storing & in het 
systeem aanbrengt en de door E veroorzaakte storing ÔW in de energie- 
integraal uitrekent. Als 8W < O0 is de toestand instabiel, Dit geldt al= 
leen voor systemen met potentiële energie (dus v = 0) 

Bernstein, Frieman, Kruskal ern Kulsrud, Proc. Roy. Soc. A 244, 17 
(1958) vonden dat een noodzakelijke en voldoende voorwaarde voor stabili= 


teit is dat de energie variatie 


nw = 4} | dv ea ed q* E + yp (v.5)? + (V,E) (6.Vp) ll > 0 
(17) 
met 


Q = Vx{EXB), | (18) 


vooropgesteld dat E‚n = 0 op de rand van het volume V , n is de nor- 
maal vector. Het zij hier opgemerkt dat (17) een tweede orde integraal 
is. Dit is zeer belangrijk omdat een eerste orde integraal voor kracht= 


vrije velden altijd &WÌ = 0 geeft; immers 








B.B B.VXEXE 
Sul = | =d e= | av . (19) 
4m 4m 
v v 
Met Velaxbl= b, Vxa'-Vxb . a volgt uit (19) 
(VxB) . (£B) Vsl (EB) xB |] 
sul = L IV + | dv . (20) 
4m 
v Ni v 
Met de stelling van Stokes wordt de laatste integraal: 
Ì ee 5 
Aivstepssl.a-gfepsr.e.0, (21) 


Vv ö S 


oee 


omdat &£ = 0 op S . Er geldt dan voor sui ; 


ôw! „lernen. (22) 
Vv 


Voor een krachtvrij veld geldt (12) en dus is &WÌ = 0 , Het toont aan 
dat krachtvrije velden met een minimale energie-inhoud zijn, maar zegt 
niets over de stabiliteit. Daarvoor moeten de tweede prde termen worden 
meegenomen. Een lineaire stabiliteitsanalyse kan in dit geval niet worden 
toegepast. 
(Zie ook: Voslamber en Callebaut, Phys. Rev. 1962, vol. 128, 2016.) 

Als v #0 ‚ kon door Frieman en Rotenberg (Rev. Mod, Phys. 1960, 
32, 898) worden aangetoond dat als E de uitwijking uit de evenwichts- 
toestand is en F(E) de kracht die t‚g.v. de uitwijking 5 op het sys- 


teem werkt een voldoende voorwaarde voor stabiliteit is dat 
me |Een so. (23) 


In een krachtvrij cylindrisch filament zorgt het axiale veld B, er 
voor dat de stroom niet t.g.v. het B, veld kan pinchen. Deze toestand 


P à 5 En : Ï 
is stabiel tegen zuiver radieïe storingen. Het veld B, zo — en B, « Ee 


6 ® 
Als we nu de straal vergroten met een factor a dan neemt B, ster= 
ker af dan Bb en zal het pinch-effect er voor zorgen dat de buis weer in 
de evenwichtsstand komt. Hetzelfde geldt voor een kleine compressie, Het 


axiale veld werkt dus stabiliserend tegen pinches, 








6.2. Annihilatie van magnetische energie en reconnectie van veldlijnen 
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Bij het optreden van magnetohydrodynamische instabiliteiten kunnen 
hoge stroomdichtheden voorkomen, waardoor de dissipatie zeer sterk toe- 
neemt. We zullen daarom eerst de annihilatie van magnetische energie als 
gevolg van het eindige geleidingsvermogen van het plasma bespreken en 
later op de eigenlijke instabiliteit terugkomen, Beschouwen we een neu- 
trale laag met aan weerszijden magnetisch veld, dan is de stroomdichtheid 


j gegeven door 


le DAE 
De diffusiesnelheid Vp is gegeven door 
Vn = (2) 


D Lm ok 
waarbij 4 de “dikte' van de neutrale laag is. 
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figuur 8 


De karakteristieke tijd waarmee de diffusie plaatsvindt is 


Tp = 4m 0 ge , (3) 


Stel voor de eenvoud dat de stroomdichtheid een Gaussfunctie is 


2 


p zg? 
j * Ji ®@ N (4) 


dan is de energiedissipatie per eerheid lengte langs de “sheet” 


2 
+0 jo? °T JE de ÍE 
2 0 2 0 
== d EE an ra} £ TE me /E 
ve fba fe Nae dt. (5) 
_—_00 00 


ieden. aldhone eden ndi 


Met 


; „AE ko (6) 
ee mel SS 
bn | dx wel 4m £ 
wordt (5): 
en 
Bo? Je Bo ° " 
We mer 2% ie 
Bo° 
Dus de magnetische energie ri die met snelheid Vo zowel van =® als 


ven +e (factor 2) naar de neutrale laag diffundeert is precies de hoe= 
veelheid gedissipeerde energie, Als S de effectieve oppervlakte van de 


stroomlaag is geldt voor de energiedissipatie 


Bo ° 
Woe S Erie (8) 
Stel dat de effectieve oppervlakte S= 101° cm?, (misschien zelfs te 
groot), dan geldt met We = 10°° voor een grote zonnevlam 
Ba wm, & 10, (9) 
D 
Met Bo = 100 gauss volgt hieruit Vp 5 107 em/sec . Dit is vergelijkbaar 


met de Alfvên snelheid: 
Bo 


1e sk 3 . 107 em/sec (o= 10 Ì* g/em?) , 
„/ 4T p 





hetgeen betekent dat door diffusie in een vlakke laag net kan worden vol=- 
daan aan de energie-eis, Anderzijds moeten we ons wel bedenken dat een 
diffusiesnelheid N= 107 em/sec in een plasma met een geleidingsvermo- 


gen o= 10 ° een dikte geeft: 


En Ì de =2 
LL = ers 7 10 * em. (10) 


Het lijkt erg onwaarschijnlijk dat zulke dunne lagen kunnen ontstaan. 
De enige manier om een dikkere stroomlaag te krijgen is: 

1) OJ kleiner 5 en 

2) Vo kleiner. 

We zullen in de volgende paragraaf zien dat er oorzaken zijn waarom O 
kleiner is. vj kleiner maken kan echter alleen als Bo groter is. 
Bijv. als Bo = 10° gauss kan v, = 10° cm/sec zijn. 


D 
Met 5 = 10 So = 10 wordt L= 10° em. 


In bovenstaand sheet model wordt steeds meer materie naar het 
centrun toegebracht. Als men naar een stationaire oplossing zoekt blijkt 
dat in een ééndimensionaal model niet mogelijk. In twee dimensies wel, 
De toevoer van veld en materie gebeurt dan in de ene richting en de af- 
voer in de andere. De totale hoeveelheid energie die ter beschikking 
staat is nooit meer dan in vj per oppervlakte eenheid, maar nu komt 


aiet alle energie in warmte terecht maar ook cen gedeelte in kinetische 


energie. 
y 
hen je B 
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figuur 9 


Laat weer # de dikte van de stroomlaag zijn en L de lengte. 


Behoud van massa geeft voor een incompressibel medium 


V) L = vi 4. (1) 


Als het medium compressibel is kan mi * no zijn dan geldt 


Vp to L = vj ar LL. (12) 


De snelheid vj waarmee de materie uit de laag wordt weggeschoten wordt 
bepaald door de druk in de laag. Deze is ongeveer gelijk aan de magneti- 
sche druk omdat druk-evenwicht geldt in de y-richting. 


2 


Bo 
Pi 9 Da tege (13) 
Bi Ro? 
Als in het plasma buiten de laag iT >> po geldt pa Er . 


Dit betekent dat de uitstroomsnelheid ongeveer gelijk is aan de Alfvén 


snelheid. 
Bo 
VF (14) 
Â wv 4m U ni 


waarbij u de massa van een proton is. 


Uit (11) t/m (14) volgt dan 


no L Ì 
4e are | (15) 
A 


ni ár 0 v 


De maximale snelheid waarmee dit proces ken plaatsvinden wordt natuurlijk 


weer bepaald door de diffusietijd (vgl. 3) die met (15) wordt: 


no L 
En (16) 
De 
D ni Va 
1 
Bij een compressie en 100 en L e= 10° em en V‚ mi 10° em/sec 
o A 
is T, = 10 sec , hetgeen zeker niet erg klein is. De diffusiesnelheid is 
Ad 
echter vergeleken met de Alfvêén-snelkeid wel klein: 
Di , Î 
ee el Ì ne ee (17) 
D 4m a 2 no 4n C L 
KAN 
Als Mo = keg het magnetische Mach getal van de instroming is en 
Â 
R, = är JL Vv, het magnetische Reynoldsgetal, geldt 
nr jé er | 
Mo = _— Ô | 
to | (R) (18) 


In het algemeen is Re groot, bijv, KR. = 10Ï! als L= 10°, u 10° 
en 0 = 10 * . Dan is Mo #10 * — 10 5 , Dus de instroming geschicdt 
zeer langzaam en is geheel bepaald door L en de Alfvénsnelheid Va 
Merk op dat ook in dit proces f% heel klein is. Vit (15) volgt 


1 
== 10 cm. 


Petschek wecs er op dat bij het proces van de annihilatie van veld- 
lijnen in het algemeen reconnectie zal optreden en dientengevolge de an=- 
nihilatie niet alleen door diffusie bepaald is, maar ook door hydromag- 
netische effecten. Dit is vooral belangrijk als de veldlijnen een hoek 


met elkaar maken, 


\ 





“diffusie” “mhd golf" 
figuur 10 


hd 14, 


In bovenstaande figuur zal aanvankelijk de diffusie er voor zorgen dat de 
scherpe hoek tussen de krachtlijnen wordt afgestompt. Later zal de sto- 
ring (en dus de afvlakking) zich als een magnetohydrodynanische golf 
voortplanten. 

In feite planten zich twee m.h.d.-schokgolven langs de veldlijnen 
voort. Deze schokgolven kunnen alleen op hum plaats blijven als het 
plasma met de veldlijnen naar de y-as wordt toegebracht. 

Petschei rekende de maximale annihilatie uit in een stationaire 
situatie, Daartoe verdeelde hij het gebied in drie delen: 

Ì. een gebied waar het plasma langs de x-as naer binnen stroomt; 

2, een gebied waar het plasma langs de y-as naar buiten stroomt; 

3. een diffusiocgebied Ky zodanig, dat voor [xl] <x en Iyl < y de 
stroming door diffusie is bepaald. 

Gebieden 1 en 2 zijn door een stationaire schokgolf gescheiden 


schok 
/ 





figuur 11 


vene - 


15. 


Als we aannemen dat in gebied 1 de snelheid alleen een component in de 

x=richting heeft, kan uit de schokrelaties voor m‚h.d.-schokken worden 
lei = = 5 — NN ni 

afgeleid dat Bz B0) met B, Bi, ‚ Dus de x=component van het 

magnetische veld blijft ongewijzigd terwijl de y=component verdwijnt. 
De schokgolf plant zich langs de veldlijnen met de snelheid 


B 
ii Se 


4m o 


voort. Dat betekent dat alleen een stationaire oplossing 





; 
mogelijk. is als het plasma in gebied Ì met een snelheid vi instroomt , 


B 








welke snelheid gelijk is aan vi = E_… , Definiëren we het tach getal 
Jam p 
als 
Vi 
Teng (19) 
v 


dan geldt vj = M Ve 
re 


Uit de continuïteitsvergelijking volgt: 


Vi WY == V2 X. (20) 


Uit het impulsbehoud in de y-richting volgt in het geval van een 


stationaire stroming: 


a 2 ee Xx 
zp P va“ xXx) tr n (21) 


De term rechts in (2í) is de Lorentzkracht CPB), geïntegreerd over de 


schok 





ï \ d B 
GD, zj) Erve, an 3 Els pe XT 
B E en 
== handen n == eten aensnd (J 
LT | d be 4n we 
omdat B over de schok constant is. 
Uit (19), (20) en (21) volgt 
ERD 
M me Bm Ed 
Ene Bn (23) 
Uit vgl. (23) volgt 
x= M |yl (24) 


Dit betekent dat de dikte van de grenslaag tussen de schokken lineair 
toeneemt met y . Hierbij is verondersteld dat vi ” vily) en o*o(y). 


Uit (26) en (24) volgt direct met (19): 


a . (25) 


Dus de materie stroomt met de Alfvén snelheid het gebied 2 uit, 
In de buurt van het x-type neutrale punt overheerst de diffusie, 


Hier moeten we weer de stromingssnelheid vj gelijkstellen san de 


diffusiesnelheid 
us, (26) 
Lm 0x 
of met (19) 
Went Been (27) 


* 
4m J Vv, Xx 
A 


ke) 


Woerdt op grote afstand het Mach=-getal bepaald door (19), dus M= E 5 


vlak bij het neutrale punt wordt de dikte van het diffusiegebied y 


bepaald door (27). Als M gegeven en v, gegeven is x bepaald door 


AN 
% Ì 
EEN 2 
E án 0 v‚ dM . (28) 
Substitueren we deze waarde voor x in vergelijking zy, dan vinden we 
% 
voor y 
* Ï 
De en (29) 


8&n 0 v, Mm? 
Á 


Voor y < y hebben we een door de diffusie bepaalde oplossing en voor 
y> y een door de golf bepaalde oplossing. Het voornaamste effect van 
dit door Petschek ingevoerde model is dat het de lengte van het diffusie- 
gebied reduceert in vergelijking met het Sweet-Parker model, Hier was de 
maximale annihilatie bepaald door het Reynolds getal ee bn 5 L Va 
In het Petschek model stelt de stroming zich ze in dat bij gegeven in- 
stroomsnelheid y bepaald is, 

Is vi groter dan het maximum, gegeven in vgl. (18), dan zal y <L 
worden en voor y > yv het proces zonder diffusie verlopen. 

Bedenk dat hoewel diffusie natuurlijk wel een rol speelt, het proces niet 
alleen door de diffusie wordt bepaald. Petschek rekent uit dat de maxima- 
le annihitatiesnelheid wordt bepaald door de verhouding L / y met L 


weer de lengte in de y=richting van de schok, 


M = pete Ne, " (30) 
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% . 
Met y gegeven door (25) vinden we voor Max 
Max 7 : Gla) 
a - ” 2 
4 Am (8m 0 VA L Wax ) 
of ook 
a evennoeen mwn (31b) 
ap 4 An (2R M 2 
m max 


Het is duidelijk dat het geleidingsvermogen wel cen grens stelt aan Max 
maar vergelijken we (31b) met (18) dan zien we dat in formule (316) Max 
slechts logarithmisch van het magnetische Reynoldsgetal afhangt en in 

(18) omgekeerd evenredig met de wortel uit B is, Hier is dus een grote 
winst geboekt met betrekking tot de annihilatiesnelbeid, Vergelijking 31b 
is een impiiciete uitdrukking voor Max . Voor grote Reynoldsgetallen 

(grote 0) is He 
In het geval van een compressibel plasma gaat de beschouwing op de- 


>> Mo « 


zelfde manier. Men vindt dan: 





M = if == (32) 
Ct tz) Am CZ Ra MN zÔ ) 


a 


Stel R = 10ÌÌ dan &n (22 M °°) 10 en we zien dat M __ “Ol, 
m m max max 
een resultaat dat later door Priest en door Parker ook wordt vermeld, 


Dit is echter veel sneller dan in het sweet mechanisme waar 


_ Rn cc 8 
Mo > R Ô 10 7? , De annikhilatie=snelhcië is een factor 30 000 
1 


groter. 
Aan het Petschek-model is uitvoerig gewerkt door Yeh en Axford (1970), 


Priest (1972, 1975), Sonnerup (1969), Parker (1973). 
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De stroomdichtheden in neutrale lagen kunnen zo groot worden dat de 
electronendriftsnelheden, Vp» vergelijkbaar worden met de ionen-, respectie=- 
velijk electronen thermische enelheden, Als dat gebeurt kunnen gemakkelijk p 
plasmagolven worden opgewekt die 
a) makkelijk hun energie kunnen overdragen aan snellere deeltjes (versnel- 

ling); 
b) het geleidingsvermogen sterk reduceren waardoor grotere dissipatie op= 
treedt (verhitting). 


Het opwekken van plasmagolven gebeurt zodra 


' (Ii) 


> 
Yp Yphase 


waarbij v de fasesnelheid van plasmagolven is. In een magnetisch plas=- 


ma kunnen Ee verschillende golfmodi worden aangeslagen; we zullen ons 
daarom beperken tot een niet magnetisch plasma. 

Hierin kunnen zich twee soorten longitudinale plasmagolven voortplanten, nl. 
1. De Langmuir golven en 

2. De ion-acoustische golven. 


Voor Langmuir golven is de dispersievergelijking 
: 





2 2 sar | d Ve E 
w (kk) = We + 3 Vre k | ag Woe + nT (2) 
pe 
bre? n_ 1} 3 a 
waarin w = ( = | = 5,65 Xx 10* n ° sec Ì de electron plasma- 
pe \ a e 
Bel k Te : 
frequentie is, bt ad. ier de electronen thermische snelheid en k 
e 
3 Vr & ië 2, 
het golfgetal. In het algemeen is w D> ze zodat w (k) Zw . 
pe 2 Ü oe pe 
De fase snelheid van de Langmuir golven 
“pe | 

Ne tn B (3) 

Voor ion-acoustische golven is de dispersievergelijking 
n Vs k 3 Vri k? 
ener rr ant (4) 
VI+ d)? Pi 
/k T Vre 
waarin Nn Ek ionengeluidssnelheid is, See Od Debye 
i i 


lengte, de ionenthermische snelheid en de de ionenplasmafrequentie, 


Vn: 
Ti 


ON 
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We gaan niet verder op het gedrag van deze golven in maar vermelden hicr al- 
leen dat ion-acoustische golven zich alleen maar kunnen voortplanten als 

Te >> Ti . Beide typen golven kunnen zeer sterk gedempt worden door de 
interactie met deeltjes die ongeveer dezelfde snelheid hebben als de golf 
(Landaudemping). Anderzijds kunnen deze golven natuurlijk net zo gemakkelijk 
worden opgewekt door deeltjes die ongeveer even snel zijn als de golf. 
(Cerenkov conditie). Demping of emissie hangt af van de gradiënt van de 


verdelingsfunctie f(v) bij de fasesnelheid 


df | > 0 > emissie , 
v 


Bv 
ph 
(5) 
zE | < 0 + demping . 
v 
ph 


In onderstaande figuur is het verband tussen w en k voor beide golftypen 


aangegeven. 


Langmuir 
pe 


pi 


Ion acoustisch 


figuur 12 


Als de stroomdichtheid toeneemt zal eerst een moment optreden dat 
Ty Ty zodat ion-acoustische golven worden opgewekt, 
Het belangrijkste effect hiervan is dat het geleidingsvermogen sterk 


gereduceerd is 


2 
ü) Vv, 
og = PV Z 10 w ie ” (6) 
TT Vr pe Vp 
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x Te 
s 10 * 


De vrijgekomen energie is dan vele malen groter dan zonder plasmaturbulen- 


Dus o° > 100w omdat v 
pe 


tie (d.i, een toestand zonder een spectrum van plasmagolven). 
Laten we nu eens nagaan hoe groot de stroomdichtheid moet zijn opdat 


anomale geleiding een rol kan spelen en of deze toestand kan voorkomen in de 


zon. 
. . _ ec dB _ cB mt a E 
Als j ne vj en j = Zr 3 ) betekent Vj Te dat 
sn sl. 0) 


An ne £ v 
s 


Als in een stroomlaag drukevenwicht heerst en B(x=0) = 0 geldt 


B? 
Ep” ZnkT. (8) 


Uit (7) en (8) volgt 


Vv 
Bl = 3,5 X 10° , (9) 
Met T = 2x 10° K en ie 1,5 X 107 em/sec volgt uit (9) 
BA = 5 Xx 10°, (10) 


Dus als B = 100 gauss geldt #4 = 50 cm, Er zijn sterke veldgradiënten 
nodig om ion-acoustische turbulentie te krijgen. De reden hiervoor is dat de 
dichtheid in de neutrale laag bepaald is door de magnetische druk, Eigenlijk 
zou de dichtheid veel lager moeten kunnen zijn. Dit is in een neutrale laag 
niet mogelijk, maar wel in een stroomlaag waar nergens B = 0 en de dicht- 


heid vrijwel onafhankelijk is van de veldsterkte, 


Het electrische veld in een stroomlaag versnelt een deeltje, maar als 


E klein is verliest het versnelde electron z'n juist gewonnen impuls door 


2 


botsingen. De wrijvingskracht t.g.v. botsingen is evenredig met v 2 als 


Vv > v Dit betekent dat er een zeker aantal deeltjes is dat steeds 


Te * 
meer impuls wint en dus “runaway” wordt. 
Men spreekt van het runaway of Dreicer electrische veld als het veld 


waarvoor de deeltjes als het ware wegvliegen (runaway) 


E == Be GI) 


Bedenk dat als E < E) er altijd een aantal deeltjes zullen zijn die snel 


genoeg vliegen om versneld te worden. 
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Deze snelle deeltjes vormen een bundel met v > Vre CP zullen dus 
Langmuir golven opwekken. Al veel eerder kunnen ion-acoustische golven wor- 


den opgewekt. Dus het geleidingsvermogen is al gereduceerd en dus 


4 


E = >> zE, (12) 


a fe. 


on 
& 
o 


als j bepaald is door uitwendige condities zoals het magnetische veld. 
Bovendien veroorzaakt de ion-acoustische turbulentie een verhitting 
waardoor E‚ daalt en de verhouding E / E) zeer sterk stijgt. Het resul- 


D 
taat is een bundel snelle electronen die Langmuir golven opwekken bij een 


golfgetal kg 2E, waarbij v de snelheid van de electronen is. 

Zolang de energiedichtheid van de plasmagolven laag is zal door ver- 
strooiing van de plasmonen aan de Debye wolk rond de ionen de energie we 
van de plasmonen gelcidelijk naar lagere golfgetallen (grotere golflengten) 
verschuiven, 

Dit proces gaat net zo lang door tet een kleinste golfgetal ko wordt 
bereikt, waarbij de plasma-interacties zodanig gecompliceerd worden dat ver- 
dere verschuiving niet meer mogelijk is. In principe zou ko = ze kunnen 
zijn, maar reeds voor die tijd treden allerlei verzadiginsprocessen op. 

Als reeds ver voor ko wordt bereikt de golfenergie wordt uitgestraald 
(plasmastraling) spreckt men van stralingsturbulentie, Is het plasma een in- 
efficiënte straler, dan zal de bij Ee ingevoerde energie zich bij ko op= 
hopen totdat Wi (ko) zo groot wordt dat men niet meer van zwakke plasma- 
turbulentie kan spreken. 

Als 


9, 
We (ko) 


Ed 


p 

mann A. (13) 
nk T 

spreekt men van sterke plasmaturbulentie. Het energietransport vindt dan 

plaats van kleine k naar grote k , Voor een bespreking van dit effoct 


verwijzen we naar Tsytovich, 1970 en Kaplan & Tsytovich, 1974, 


plasmonophoping 
| (sterke turbulentie) 


te 


figuur 13 
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De figuur op blz. 21 geeft het encergietransport weer, het energie- 
transport van K, naar ko “ontmoeten! de plasmonen steeds minder deeltjes 
met steeds hogere snelheid waarmee ze in interactie kunnen komen, De plasmo- 
nen zullen dus een beetje gedempt worden (Landau demping), waardoor de snel 
le deeltjes nog sneller zullen gaan bewegen. Er treedt dus versnelling op. 
Bedenk dat 


u. 
v' = ei > RS, = Vv (14) 


als k < kg . 


De conclusie is dat bij zwakke plasmaturbulentie er een klein aantal 
deeltjes is dat tot zeer hoge snelheden kan worden versneld. 
De overgrote meerderheid van de deeltjes (thermische plasma) bevindt 


zich rond wv of in de k-ruimte bij k 7 aad Ke . Zodra het 


Te De Vre 
moment optreedt waarop de plasmaturbulentie zo sterk wordt dan van energie= 
omkering kan worden gesproken, ‘stroont' energie terug van de energierijke 
plasmonen rond ko naar de plasmonen rond Be . 

Als we mn k £ kan het thermische plasma worden verhit. Men spreekt 
dan van verhittingsturbulentie door plasmonen (Langmuir plasmonen). Er 
treedt dan evenwicht in tussen thermische plasma en plasmaturbulentie en een 
thermalisatie zorgt er voor dat geen verdere deeltjesversnelling meer op- 
treedt. 

Onderstaande figuur laa: het energietransport zien voor sterke plasma- 


turbulentie. 


W evenwichtsspectrum 





figuur 14 


Er stelt zich dan een evenwichtsspektrum in volgens in volgens Kolmogoroff, 


Voor plasmaturbulentie'is dit spektrum gegeven door 


Ed 
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ue iS (15) 


(Een Kolmogoroff spectrum) 


Als conclusie kuanen we stellen dat bij voldoend hoge stroomdichtheden 
in het vlamplasma achtereenvolgens efficiënte versnelling van deeltjes plaats 
vindt gevolgd door verhitting van het plasma zoals dit ook is waargenomen ín 
zonnevlammen, 

Men vermoedt dat de stralingsprocessen pas belangrijk worden als de 


sterke verhitting van het plasma plaats vindt. 


